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1 Introduction 



Abstract 

A halo model with heavy relic neutrinos N belonging to a fourth 
generation and their annihilations in galactic halo may explain the recent 
^ ' evidence of diffused gamma (GeV) radiation around galactic plane. We 

P^ , considered a neutrino mass in the narrow range {Mz/2 < tun < Mz) 

t^ ' and two main processes as source of gamma rays. A first one is ICS of 

ultrarelativistic electron pair on IR and optical galactic photons and a 
second due to prompt gammas by vr" decay, leading to a gamma flux 
(10"'' - 10"'' cm-'^ s"i sr-^) comparable to EGRET detection. Our 
^^' \ predictions are also compatible with the narrow window of neutrino mass 

■"■^ — ■ 45 GeV < TUN < 60 GeV , required to explain the recent underground 

1-^ ' DAMA positive signals. 

O 

^ ' The recent observations of EGRET telescope show a diffuse 7 - ray emission 

in the halo of our galaxy [1]. Actually models trying to explain a gamma halo 
range between 

a) a high galactic latitude distribution of high energy cosmic ray sources (fast 
\^ ' running pulsar " Geminga" like) [1] [2] , 

b) collisions of high energy protons (tens of GeV) in molecular clouds (mostly 
H2), whose formation is favored in galactic halos [3], 

c) a cold dark matter scenario with neutralino, the lightest supersymmetric 
particle in the minimal SUSY model, whose annihilations in heavy fermions 
{cc, bb, tt) and bosons {W^W^ , ZZ, gg) could lead to gamma secondaries 
emission [4], 

d) Inverse Compton Scattering (ICS) of interstellar photons off cosmic ray 
electrons with a spectra harder than previous predictions [5] . 
A heavy Dirac neutrino of a fourth generation in a Cold Dark Matter (CDM) 
model, at masses mN > Mz/2, may offer an elegant solution to the gamma 
ray signal detected, though it can not solve all the dark matter problem in the 



halo. 

Neutrino annihilations at high galactic latitude could produce 

1) 7 rays by ICS of rclativistic electrons (primaries or as secondary decay prod- 
ucts of heavier leptons in the annihilation chains NN -^ l^l^ -^ e+e^) onto 
thermal photons {IR, optical) near and above the galactic plane, 

2) direct gammas by neutral secondary pions decay. 

The estimated flux we derived here is roughly close to EGRET results in the 
hypothesis of a smooth and homogeneous galactic halo. 



2 The heavy neutrino model 

LEP I has fixed sever constrains on the number of "light" {m^ <C Mz/2) neu- 
trino families from Z width data. However there is no experimental prohibition 
on the existence of an additional heavy neutral lepton with mass Mn > Mz/2. 
Such a heavy stable neutrino belonging to a fourth fcrmion family was intro- 
duced nearly 20 years ago as a CDM candidate [6], [7]. An apparently simple 
fourth generation model is not so easy to build. Anyway models predicting a 
fourth heavy stable neutrino can be found in [8], [9], [10]. In an expanding uni- 
verse scenario heavy neutrinos decouple from a thermal equilibrium condition 
when global temperature drops below their rest mass energy (T < iiin) and 
weak interactions become too slow to keep neutrinos in equilibrium with the 
cosmological fluid. From this moment the only change in neutrino density is 
due to cosmic expansion. The relic abundance is given by [11] 
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where g* = Ntos + |-^^/erm is the number of effective degrees of freedom 
at temperature T, gs is the number of particle spin states Mp is the proton 
mass, {(y(3)f is the thermally averaged annihilation cross section at freeze out, 
n-y = 0.24T^ is the cosmic photon number density. 

Annihilations of heavy neutrinos in the universe happen through two main 
channels 

channel 1) NN -^ ff if Mz/2 < tun < mw 

where the cross section decreases as Arnj^^ / {Arnj^ — -^i)^ for growing tun 
channel 2) NN ^ W+W^ if niN > mw 
with a cross section growing like mf^ for increasing mj\j [12]. 
As pn (X a~^, neutrino relic density exhibits a maximum (with Pmax/Pc ~ 



lO^^/i^^) near mpf r^ Mw and then starts to decrease as m^ in the mass 
range uin > ttiwi without reaching the critical value il = 1, at least in the 
mass range where Standard Model may be applied {mN < 1 TeV) . 
Clustering during galactic structure formation determine an increase in neu- 
trino density [13], [14], [15] that in the central part of the galaxy could be as 
large as 5 -^ 7 orders of magnitude (the exact value depends on the other CDM, 
HDM densities and masses) . 

A spherical halo around our galaxy made of heavy neutrinos is the model pro- 
posed in order to explain 7 emission observed at high galactic latitude. In 
a galactic halo, neutrinos with a higher density distribution could annihilate 
again leading to a flux of ordinary particles beyond the galactic plane, poten- 
tial sources of high energy radiation. 

Constrains on neutrino mass come either from cosmological data (not too high 
pollution of e+e~ cosmic rays is observed in the range Mz < niN < 300 GeF 
[11]) or DAMA detector, where recent signals could be attributed to heavy 
neutrino in the mass window 45GeV < niM < 50Gey[15]. 

3 Neutrino annihilation products as gamma ray source 

3. 1 Relativistic electron pairs: ICS on the galactic interstellar radiation field(ISRF) 

Heavy neutrinos could directly annihilate in relativistic electron pairs (either 
prompt ones or born through secondary decay processes of heavier particles 
/i, r, as in channel 1 way) . Channel 2 leads to electron pairs through leptonic 
decay of W {NN -^ W~^W~ -^ l~^l~). Electron pairs may be generated even 
by W, Z hadronic decay through charged pious and neutrons production. 
Electrons and positrons are trapped by galactic magnetic field, and propagat- 
ing through the Galaxy loose either "memory" of their "place of birth" as 
well as energy for bremsstrahlung, synchrotron or ICS. These processes de- 
termine a broadening of different electron "lines" {NN -^ l^l^ -^ e+e^), so 
that e~ (e"*") spectra (even considering electrons and positrons that come from 
hadron decays) are at final stages described by the consequent approximated 
power law J = KE^"' (where K is a normalization constant). Numerical sim- 
ulation of NN annihilation performed with the package PYTHIA 5.7 [16] with 
suitable modifications to include a fourth generation of fcrmions, show that 
such NN relic electron fluxes are considerably lower than observed neighbor 
galactic background (in the range of masses AhGeV < rriN < Mz). Such a 
cosmic ray input can not be used to confirm or refute heavy neutrino presence 
in galactic halo. 

ICS of "soft" background photons (isotropic CBR or anisotropic infrared and 
optical interstellar radiation field) off relativistic electrons created out of the 



galactic plane is a possible source of radiation. Energies order of magnitude is 
fixed by the ICS characteristic relation E^ = A/'Seph{Ee/meC^)'^ , where eph is 
the target photon energy. 

We excluded here collisions with microwave photons which would require too 
large neutrino masses {rriN > ITeV). No clear theory is still available for 
such a heavy particle. ICS on IR and optical photons needs respectively 
Ee > 50GeV and Eg > lOGeV, and is more efficient in gamma ray pro- 
duction. The Galaxy is a disk-like radiative source of radius ^ 15 kpc, so the 
interstellar radiation field has a vertical extent of several kpc. We assumed 
this radiation to be represented by the obvious scaling law 

l + r-'/ai^ 

where r is the distance from the galactic plane, and a-y ~ 10 kpc is the 
characteristic length of interstellar radiation distribution in the Galaxy. Pho- 
ton density could be considered roughly constant in a region of radius a^ [15]. 
An electron distribution {KE^"") interacting by ICS with photons at energy 
eph and density Uph generate radiation whose intensity is [17] 

ME,) = lKa,np,ar (^) '" '''' E-i^+^^Z^ (3) 

where Hph is the target photon background density, eph its average energy 
and (Tt is the Thomson cross section. 

Gamma intensity has been calculated for totv ~ 45, 50, 100, 300 GeV. The 
largest flux has been obtained for ICS on optical photons. 
For rriN = 50 GeV the calculated flux is 
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and for ttin = 100 GeV one flnds 
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dSdtdndE, ^^'\GeVj \10kpc, 

(5) 
A{ip) is the adimensional integral of interstellar photon density along the 
line of sight L, defined by the angular coordinate ip (angle between L and the 
direction of the galactic centre). A{ip) is of few unities and corresponds to 



O'l J line of sight U + ^IVj /O^J 

Gamma intensity due to infrared background is less abundant than optical 
photons as a consequence of the spectral power law E^°' . 
Assuming an average NN clustering p^ /pj?"™ = 10®, the flux obtained for 
two values of niN 
1)%{E > 1 GeV) 2± 4 • 10'^ Ai^) [j^) cm-^ s'^ sr-\ {niN o^ 50 GeV), 

2)^^{E > 1 GeV) 2± 6 • 10"^ A{^P) (j^) cm~^ s^^ sr~^ (m^ ~ 100 GeV), 
is comparable with EGRET observations: 
^^{E > 1 GeV) ^ 8 • 10-^ cm-^ s' 
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Additional tests of this model could be obtained with the nearly detectable sig- 
nal of hundreds KeV radiation in the halo with flux J^ ~ 10~^ cm~'^ s^^ sr^^ 
at peak energy E^ r^ SOOkeV as well as a flux Jx — 0.3cm~^ s~^ sr^^ at 
Ex ^ ^ KeV due to ICS of tens of GeV and GeV electrons with GBR. An 
addtional parassite radio background arises at high galactic latitudes due to 
synchrotron losses of the same electrons at E^ '^ 10 GeV, with typical density 
flux 

.,„ . 5 . 10« (^) (^) -' (5^1^) "... (7) 

and average frequency 

.^,^(^-iGH^(^)(^L_y (8) 
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where we used as characteristic scale for the magnetic field B = 1 /iG, and 
as optical background energy density Urad = 0.2 eV cm^^ . 



3.2 Annihilations in gamma photons 

Gamma radiation could also be produced in neutrino annihilations due to 
neutral pions secondaries in Z, W hadronic decay. This kind of emission does 
not need to introduce any kind of radiative background and is directly related 
to neutrino distribution in the halo. 
Photon flux is described by the following expression: 



47rm^ ^ "-C/ J line of sight 

where ^ is the angle between the hne of sight and the galactic center, 
p{r) is heavy neutrino density as a function of galactocentric radius, and 
^jCTiW^j^ counts all possible final photon channels {^^) which could con- 
tribute to gamma photons emission. The integral of neutrino density along 
the line of sight L depends on the halo model chosen for dark matter distribu- 
tion, which is generally described as 

^(^) °^ (i:)7[i + (i:)«](/3-7)/a '^^^^ 

The simplest density profile is described by an isothermal sphere with 
a = 2, /3 = 0, 7 = 0) and length scale a > 10 kpc. 

In the spherical model the square density integral leads to an adimensional 
intensity /(V") 

this intensity has a characteristic behaviour which is maximum in the di- 
rection of the galactic center (i/; = 0), and then decreases ior Q < ip < tt, but 
it doesn't vary more than a factor ten with the angular coordinate. 
At high latitudes / {if)) is generally of order unity. 

Models with a singular behaviour towards the galactic center or which pos- 
tulate a clumpy distribution of dark matter could contribute to enhance the 
total gamma flux in the halo, but we shall neglect them here. 
Monte Carlo simulations of neutrino annihilations [16] have been also used to 
compare EGRET flux, showing that it is possible to extrapolate a power law 
for gamma spectrum. 
An approximated integral flux for ttiat ~ 50 GeV and a ~^ 10 kpc is roughly 

<^^>Q-IQ-'^ I{'ip)cm-'^s-^sr-^, (11) 



while for mj^ ~ 100 GeV at a '-^ 10 kpc 

$^ >4-10"'^/(V')cm"2s"^sr"\ (12) 



In conclusion there are at least two independent processes able to solve 
the puzzle of a GeV gamma halo by the role of a cosmic relic DM made of 
fourth generation heavy neutrinos. The neutrino masse are constrained into 
a narrow energy window (45 GeV < ttin < 60 GeV) , in order to combine at 
once the DAMA data and the other underground detector. This reality will 
be soon confirmed or excluded by LEP II search for e~e~^ -^ NNj events [19] 
in this energy band. 
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